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The one-dimensional spin-1/2 XXZ model in a mixed transverse and longitudinal magnetic field 
is studied. Using the specially developed version of the mean-field approximation the order-disorder 
transition induced by the magnetic field is investigated. The ground state phase diagram is ob- 
tained. The behavior of the model in low transverse field is studied on the base of conformal field 
theory. The relevance of our results to the observed phase transition in the quasi-one-dimensional 
antiferromagnet CS2C0CI4 is discussed. 

I. INTRODUCTION 

The effects induced by magnetic fields in low-dimensional magnets are subjects of intensive theoretical and experi- 
mental researchi2i2i!i£i&. One of the striking effects is the dependence of magnetic properties of quasi-one-dimensional 
antifcrromagnets with anisotropic exchange interactions on the direction of the applied magnetic field 7,8 ^. The basic 
model of such type of magnets is the anisotropic Heisenberg chain - so-called XXZ model. It is, therefore, important 
to study the dependence of the properties of the XX Z chain on the field direction. There are two studied cases of 
the field direction. First of them is the XX Z model in the uniform longitudinal magnetic field. This model is exactly 
solved by the Bethe ansatai^ and is studied in great details. In the second case the field is applied in the transverse 
direction. The XXZ model in the transverse field can not be solved exactly and various approximate methods have 
been used to its studjiii*i2iiiii4. The behavior of the XXZ model in the symmetry-breaking transverse field is es- 
sentially different from the case of the longitudinal field. In particular, the transverse field induces the perpendicular 
antiferromagnetic long-range order (LRO) and the ground state quantum phase transition takes place at some critical 
field, where the LRO and the gap in spectrum vanish. The phase transition of this type has been observed in the 
quasi-one-dimensional antiferromagnet CS2C0CI4 7 . In fact, the magnetic field can have both the longitudinal and the 
transverse components. For example, the magnetic field in recent neutron scattering experiments on CS2C0CI4 has 
been applied at an angle to the anisotropy axes. From this point, it is of a particular interest to study the ground 
state properties of the spin s = ^ XXZ chain in coexisting longitudinal H z and transverse magnetic fields H x . The 
Hamiltonian of this model is given by 

N N N 

H = / ^(SnSn-i-l + Sn^n+l + ^^n^n+l) ~ ~ J]] (1) 

n—1 n—1 n—1 

where 

, 9x(z)HbH x{z) 

h x(z) - -j {*) 

is the effective dimensionless transverse (longitudinal) magnetic field, J is the exchange constant and A is the 
anisotropy parameter, which is assumed to be A > — 1. 

It was proposed^ that low-energy properties of CS2C0CI4 in the external magnetic field is described by the Hamil- 
tonian Q with A = 0.25 and J = 0.23 mev. 

Evidently, in the case A = 1 the behavior of the system does not depend on the magnetic field direction and the 
model (JIJ reduces to the isotropic Heisenberg chain in a magnetic field h = \Jh 2 z + h x . In the limiting case A — * 00 
the model JU reduces to the antiferromagnetic Ising chain in a mixed longitudinal and transverse field. This model 
was investigated ini5il&, where it was shown that there is a critical line in the (h x , h z ) plane, where the ground state 
phase transition takes place. The critical behavior in the vicinity of this transition line belongs to the universality 
class of the two-dimensional Ising model. 

Thus, the physics of the model JIJ is very well understood in the case h x = and is fairly good for the cases 
h z = and A — > 00, but no detailed studies are available in general case. In this paper we study the model Q 
using the mean- field approximation, which is the generalization of the approach developed ini£ for the case h z = 0. 
This method allows us to determine the transition line with high accuracy. The behavior in low-Zi^ region will be 
considered using the conformal field theory method. 

The paper is organized as follows. In Sec. II we consider a qualitative physical picture of the ground state phase 
diagram based on the classical approximation. In Sec. Ill the mean-field approach is developed and study of the critical 
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properties of the model is presented. Scaling estimations of the gap and the LRO in low-Zi^ region are given in Sec. IV. 
The special case A = — 1 is studied in Sec.V. In Sec. VI we discuss our results in relation to the experimental data for 
Cs 2 CoCl 4 . 



II. THE CLASSICAL APPROACH 



In order to provide a physical picture of the phase diagram of the model Q we use the classical approximation, 
when spins are represented as three-dimensional vectors. The variational wave function corresponding to the classical 
approximation has a form of a simple direct product of single-site spin states^ 

= (1 + AiS+)(l + A 2 S+)(1 + AjS+)(l + A 2 S+) . . . \Ui . . .) (3) 

where Aiand A 2 are variational parameters. If A\ 7^ A 2 then the ground state is two-fold degenerated and another 
ground state wave function is 

1*2) = (1 + A 2 S+)(1 + + A 2 S+)(l + AiSt) . . . Illl . . .) (4) 

The form of the variational parameters A\ and A 2 minimizing the energy is different in the regions |A| < 1 and 
A > 1. For the case |A| < 1 they can be chosen asi£ 

A l = Ae^, A 2 = AT* (5) 

The ground state energy for this case calculated with $1 (or $2) is 

E Acos2(j) A(A 2 -1) 2 h x Acos<f> h z (A 2 - 1) 



N (1 + A 2 ) 2 4(1 + A 2 ) l + A 2 l + A 2 

Minimizing this energy over A and <f> one obtains 

6 = cos — , 

2 V /(A + 1) 2 -^ 2 



(6) 



_ / 1 + A + fe, 

A il + A-h z (7) 

Two-fold degenerated ground state at 7^ is characterized by a non-zero staggered magnetization along the Y 
direction, which plays the role of the LRO parameter 

(ss> - (-ir^S (8) 

For a given value of A the line of phase transition on (h x ,h z ) plane (the transition line) is determined by the 
condition cf> — and has a form 

f + (TW = 1 (9) 

This line separates the antiferromagnetic (AF) phase with the LRO from the paramagnetic (PM) phase with uniform 
magnetization. The transition line for A = 0.25 is shown on FigQ] 

The case A > 1 can be analyzed in a similar way. In this case two-fold degenerated ground state in the AF phase 
is characterized by non-zero staggered magnetizations along the X and Z axes. But the expression for the transition 
line is rather cumbersome and we do not present it here. The transition line in the classical approximation for A = 5 
is shown on Fig [3 

As it is knownii there is a remarkable, so-called 'classical' or disorder, line which lies in the AF region in the (h x , h z ) 
plane and is given by the equation: 

(10) 



(1 + A) 2 2(1 + A) 

The classical line is remarkable in a sense that the ground state on it is identical to the classical one and quantum 
fluctuations are missing. It was shown ini£ that the ground state of on this line is two-fold degenerated and the 
exact ground state wave functions have the product form J3J and ijlfl. 
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FIG. 1: The ground state phase diagram of the model Q for A = 0.25. The transition line between the antiferromagnetic 
(AF) and paramagnetic (PM) states obtained in the MFA is shown by thick solid line and that in the classical approximation 
© by dashed line. Thin solid line denotes the classical line 1101 and dotted line corresponds to separatrix line (see Sec. III). 



The ground state energy on the classical line for any even N is 

E _ 1 A 

N ~ ~2 ~ ¥ 

From Eqs.Q one can find that in the case |A| < 1 the magnetizations on the classical line are: 

K n> 2 1 + A 

(so = h 



For A > 1 the parameters Aiand Ai on the classical line are 



A + l 



and the magnetizations on two sublattices are 

Wn)l,2 





A\ a + 


1 


1^,2 


- 1 


2 A 12 


+ 1 



(11) 



A^-^ d3) 



Wn)l,2 

(S v n ) U2 = (14) 
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Thus, the classical approach shows that the ground state is different in the regions with |A| < 1 and A > 1. 
For |A| < 1 the classical ground state is given by a configuration, where the spins on odd and even sites pointing 
respectively at angles x an d ~X with respect to the XZ plane. For A > 1 in the ground state all spin vectors lie 
in the XZ plane with the spins on odd and even sites pointing respectively at angles tpi and tf2 with respect to the 
X axis. This means that besides uniform magnetizations along X and Z axes in the AF region there is as well the 
staggered magnetizations: in the Y direction for | A| < 1 and in both X and Z directions for A > These facts are 
confirmed on the classical line, where the classical approximation gives exact ground state. 

Of course, one can not expect that the classical approach gives accurate estimation of the transition line and 
correct description of the phase transition (critical exponents)^. Nevertheless, as it will be shown below the fact of 
the generation of the staggered magnetizations in the Y direction for |A| < 1 and in both the X and the Z directions 
for A > 1 is qualitative true. 



III. MEAN-FIELD APPROXIMATION 



Previously, the mean-field approximation (MFA) has been proposed to study the anisotropic Heisenberg chain in 
the transverse magnetic fieldi^i^. It has been established that the MFA works very well if the transverse field is 
sufficiently strong and it gives qualitative results for intermediate fields. For A > —0.5 the MFA allows to determine 
with high accuracy the critical transverse field at which the order-disorder transition occurs and to describe correctly 
the behavior of the system in the transition region. The MFA is based on the Jordan- Wigner transformation of 
spin- 1/2 operators to the Fermi operators with the subsequent mean-field treatment of the Fermi Hamiltonian. In the 
case of coexisting transverse and longitudinal magnetic fields it is impossible to reduce the model Hamiltonian Q to 
a local form in terms of the Fermi operators. Nevertheless, for this complicated case the MFA can be modified. In 
this section we develop the special version of the MFA, which remains the variational approach. This approach gives 
high accuracy in determining of the transition line and correctly describes the whole ground state phase diagram. 

At first we perform a rotation of the spins in the XZ plane by an angle ip: 

S* = a* cos tp + a z n simp 
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s n = -<sin<p + a£cosp 

SI = < (15) 



where a" are new spin-1/2 operators. 

The Hamiltonian is transformed to the form 



H' = sin 2ip + - (hx cos ^ - h z sin <p) < (16) 

where 

x = cos 2 <p + A sin 2 y> 

z = A cos 2 tp + sin 2 ip 

h = h z cos ip + h x sin y> (17) 

The angle </? is a variational parameter over which we will minimize the ground state energy. 
After Jordan- Wigner transformation to the Fermi operators a+, a n 

< = \-a+a n (18) 

the Hamiltonian (jlf)f) takes the form 

WV zN 1 + a: , 

i? f = — Y + ~ + Z-^( h ~ z + ~ ~ cosk ) a k a k 

+ 4 ^ sin fc(a^a+ fc + a_ fc a fc ) + z ^ a+a„a+ +1 a„ + i + i?f (19) 

We treat the Hamiltonian Hi in the MFA, which implies the decoupling of the four fcrmion term. The Fermi 
representation H' f has non-local form. But we note, that all terms in H' { contain odd number of the Fermi operators 
a n and, therefore, (H[) = in the MFA. This fact holds the MFA in the frame of variational principle. 

Thus, in the MFA the ground state energy E$ and the one-particle excitation spectrum e(k) have the form: 

Eo/N = + f + (^-z) 7l + i±^ 72 + i_^ 73 + z( 7l 2 - 7 2 + 7 2 ) (20) 



e{k) — y(u + vcosk) + w 2 sin k (21) 



where 



u = h — z + 2zji 

1 + x 
V = — 2z 72 

w = i^ + 2z 73 (22) 

Quantities 7 i, 7 2 and 7 3 are the ground state expectation values, which are determined by the self-consistent 
equations: 

_i_ f dk ( u + v cos k 

7i = (On) = / — ( 1 





72 = (<a„ + i) 



2vr V e(k) 

dk (u + v cos k) cos k 
2^ 7{k) 



o 



i dkwsm 2 k 

73 = (a n a n+1 )=- I ^^y- (23) 



n 
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The solution of the self-consistent equations l|23(l gives the minimum of the ground state energy l|20|) in a class of a 
'one-particle' wave functions at a given angle p. Thus, one should minimize the energy (|20|) with respect to the angle 
ip, solving the self-consistent equations l(23|l for each value of ip. This means that the proposed procedure remains 
variational one. 

The physical meaning of the angle ip is to show a direction of the total magnetization of the model 

S z = «)cosp = ( - -71 ) C0S</5 



2 

S x = K)sin^= Q-^Binp (24) 

Transforming the mean- field treated Fermi Hamiltonian back to the spin operators, we arrive to the well-studied 
anisotropic XY model in a longitudinal magnetic field^ 

H xy + & + ~ « X] < ( 25 ) 

The model (125(1 has a transition line defined by the equation 

u(h x ,h z ,A) = v(h w ,h z ,A) (26) 

which separates the region u < v with the LRO represented by a staggered magnetization from the region u > v, 
where there is no LRO except the uniform magnetization l(24|l . The transition line h zc (h x , A) is determined by the 
numerical solution of Eas. l|23l) . l|2t)|) with the minimization of the ground state energy over angle p. The transition 
lines in the MFA for A = 0.25 and A = 5 are shown on FigQ]and FigElby thick solid lines. 

It is well known 18 that the critical properties of the model (|25|l belongs to the universality class of the two- 
dimensional Ising model. This means that in the MFA the gap is closed near the transition line linearly with the field 
and by the 1/8 law for the staggered magnetization. 

The MFA also shows, that for |A| < 1 {w > 0) the model has a staggered magnetization along the Y axis 

\w 2 (v 2 -u 2 )] 1/8 

-ir-S* = L ) n( '{ (27) 

while for A > 1 (w < 0) the staggered magnetizations exist along the X and the Z axes 

\w 2 (v 2 ~u 2 )} 1/8 

(( _ irK) _ K^- 2 )]' /8 5i „ y (28) 

y/2{v-W) 

The validity of the 2D Ising type of the critical properties of the model (JTJ in the vicinity of the transition line has 
been checked by the density matrix renormalization group (DMRG)A£ calculations of the staggered magnetization 
and the gap. The staggered magnetization is computed as^ 



st 



jf (olEH 1 )" 5 "! 1 ) « = (*.v.«) ( 29 ) 



where |0) and |1) are two lowest energy states. These states are degenerate (at N — ► 00) in the ordered AF phase. 
Therefore, in the AF phase the gap is given by the second excited state, while in the disordered PM phase the ground 
state is non-degenerate and the first excitation determines the gap in the spectrum. We have performed the DMRG 
calculations using the infinite-size algorithm and open boundary conditions and the number of states s kept in the 
DMRG truncating procedure is up to 25. We estimated the relative error due to DMRG truncation from difference 
between the data computed with s = 25 and those with s = 20 for chain lengths N = 202. The estimated relative error 
is of the order of 10~ 5 which is sufficiently small for accurate estimates for the gap and the staggered magnetization. 
As an example, on Figs |3 and 01 we show the plots of (A/J t ) 8 and the gap m versus h x in the vicinity of the transition 
point for A = 0.25 and fixed h z = 0.67 (these parameters are related to those for the antiferromagnet CS2C0CI4). 
A good linearity of the plotted data definitely confirms the 2D Ising character of the transition line. The excellent 
agreement between the DMRG and the MFA results on Fig|3]and FigQJshows high accuracy of the MFA. For example, 
the critical field h x estimated from the DMRG results differs from that obtained in the MFA within 0.04%. 
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1.349 



FIG. 3: The staggered magnetization near the transition line as a function of h x for A = 0.25 and h z = 0.67. 



gap 




FIG. 4: The gap as a function of h x near the transition line for A = 0.25 and h z = 0.67. (Solid line is the MFA, triangles are 
DMRG results extrapolated to the thermodynamic limit). 

In addition to the transition line denned by Ea. l|26l) . the Hamiltonian i|25[l contains another special line defined by 
the equation 

u 2 + w 2 = v 2 (30) 

This line separates the so-called 'oscillatory' region u 2 + w 2 < v 2 (lying totally in the AF phase), where spin 
correlators of the model l|25[) have oscillatory behavior with an incommensurate wavelength depending on the model 
parameters (h x , h z , A), from the region without such oscillatory behavior of correlators^. The line (|30|l is nothing 
but the classical line of the model JJJ. Remarkably, the MFA gives the exact ground state on the classical line. On 
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this line the solution of Eas. H23|) has a simple form: 



• 2 pi',— , . 

Sm<p= — P T^A (31) 



where 



The values of 7$ for | A | < 1 are 



(32) 

4 1 + A ^ ; 



and for A > 1 



1 V 1 ~P P /Qg ^ 

71=2 2 — ' 72=73 = -^; (33) 

7i = o ~ o /1 ' 72 = -73 = -rr, r (34) 

2 2y/l -p 4(1 -p) 

The ground state energy is given by {TTJ. Substituting Eas.l|3*T )l -l{3l J) into Eas.pi |l .p7 jl .p8 |l one can check that the 
magnetizations on the classical line in the MFA coincide with those given by Eas.qi2 [I .Q14 [l . 
The gap on the classical line in the MFA lies at k = ir and equals 

m = 1 - I - P, |A| < 1 

1 - 2p - pA 



m = ^r^j (1 I P) 2 ■ A>1 ( 35 ) 

It is necessary to note that the elementary excitation in the AF phase can be regarded as a domain wall between 
the two AF ground states. In the cyclic chain these excitations are created in pairs in contrast to open chains. In 
Ea. (|19fl the end-chain correction term is omitted, and the spectrum Eq.JSTJ determines the gap for the open chain 20,21 . 
Therefore, in the AF region Eqs.|J2U,£3J §i ye a nau " °f the gap for a cyclic chain. 

It is worth to mention one more special line on the phase diagram, so-called 'separatrix', defined by the equation 

uv — v 2 — w 2 (36) 

This line separates the region uv > v 2 — w 2 , where the lowest excitation has momentum fc m ; n = ir from the region 
uv < v 2 — w 2 (situated entirely in the AF phase), where the lowest excitation has momentum, depending on the 
model parameters as 

uv 

COS rZ m [ n — ~ ~ 



— V 



For any A, the transition line h zc (h x , A), the classical line h zc \(h X) A) and the separatrix h zs (h x ,A) lie in the 
following sequence (see FigQ] and Fig|5$ 

h zs (h x ,A) < h zcl (h x , A) < h zc (h x ,A) (37) 
All these lines meet each other in the only point F (h x — Q,h z = 1 + A). 



A. The point F 

The point F (h x — 0, h z — 1 + A) is the special boundary point, where all special lines terminate. The ground state 
in the point F is saturated ferromagnet. Near the point F (h x <C 1) the fermion density is small and the mean-field 
treatment of the four fermion term in Ea. (|19fl gives the accuracy, at least, up to ft 4 . We omit here intermediate 
calculations and give the final expressions for the special lines. The transition and the separatrix lines near the point 
F have the form: 

h zc (h x ,A) = h zcl (h x ,A)+m cl + 0(hl) (38) 
h zs (h x ,A) = h zcl (h X: A)-m cl + 0(h 6 x ) (39) 
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where the behavior of the classical line h zc \(h x , A) is given by Ea. (|10|) 

h 2 h 4 

h zcl (h x ,A) = l + A 



4 32 (1 + A) 

and m c \ is the gap near the point F on the classical line (see Eq.J^Hl)) 



O {hi) 



1- A 
128 V 1 + A 



0(hl 



(40) 



(41) 



As one can see the difference between three special lines near the point F is very small, of the order of h 4 ,. 
The expressions for the gap are different to the left and to right of the separatrix line: 

1 - A 



4%/2 



1 + A 



\/h zcl {h x ,A) - h z , h z < h z 



m = \h z - h zc {h x , A)| , 



h z > h z 



(42) 



The linear behavior of the gap in the vicinity of the transition line confirms the 2D Ising universality class of the 
transition line. 

The staggered magnetizations in the vicinity of the point F vanish on two lines: on the transition line and on the 
line h x = 0. The Eqs.(|ll|), d3> (EHJ near the point F reduce to 



(S*) = {-l) n B 



(43) 



2 



1- A 



1 + A 



1/4 



1/8 



for |A| < 1 and 



(S x n ) 



4 v ' 



i-iy 



KB 



(44) 



for A > 1. 

To validate our analysis in the vicinity of the point F one should also estimate the effect of the part of the 
Hamiltonian H' in Ea. (|16|l . omitted in the MFA. Near the point F the angle ip ~ h x /2 and, therefore, these terms in 
H' are small and can be taking into account as perturbations. The corresponding perturbation theory contains only 
even orders. The estimate of the second order shows that the contribution of these terms to the ground state energy 
and to the gap is of the order of hi and h x (h zc — h z ). This accuracy is sufficient to confirm the above equations. 

We note that in the limit A — > oo the point F transforms to the so-called multicritical point with macroscopic 
degeneracy of the ground stated. 



B. The point A 

In the case h z = the model reduces to the anisotropic Heisenberg chain in the transverse magnetic field, 
which was studied inAiii&ii. At some value of magnetic field h x A(A) this model undergoes the transition from the 
antiferromagnetic state to the paramagnetic gapful state. We denote this transition point by the point A (see FigsQ] 
andEJ. 

To study the behavior of the system in the vicinity of the point A we shall follow the arguments ofi&, where the 
point A was analyzed in detail for the special case A — > oo. For h x = /i x a(A) and small longitudinal magnetic field 
h z we rewrite the Hamiltonian in the form 

H = H Q + V 

Ho = ^y^X^n^n+l + + ^^n^n+l) ~ ^xa(A) ^ 

V = -h z J2s z n (45) 
The unperturbed Hamiltonian Hq describes the transition point A, where the spectrum is gapless 13 . 
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As was mentioned above, the MFA does not give the transition point ]i x a(A) exactly (though with high accuracy), 
but it correctly describes the character of the transition. Therefore, we can use the MFA to determine how the small 
perturbation V opens the gap. [We note, that the minimum of the MFA ground state energy at the point A is achieved 
at the angle <p = ir/2 and two 'bad' terms in H' in Eg. (|16f) disappear at this point.] 

To study the perturbation theory with small perturbation V one can repeat step by step the analysis done for 
the special case A — > oo of the model H45(! 16 . As a result one finds that the perturbation theory contains infrared 
divergencies, which are absorbed in the scaling parameter y = h 2 z N . The perturbation series for the mass gap has the 
form 

m = a(A)h 2 z + g(A,y)h 2 z (46) 

where a(A) is the second order correction and g (A, y) is the scaling function accumulated all high-order divergent 
terms. The scaling function g(A,y) in the thermodynamic limit (y — > oo) tends to some finite limit. So, the gap 
takes the form 

m=[a(A)+ 5 (A,oo)]/ l 2 (47) 

From the last equation we see that the gap is proportional to h 2 , but the factor at h 2 is given not only by the 
second order correction a(A) but by all collected divergent orders of the perturbation series. 

For a fixed value of A the behavior of the transition line in (h x , h z ) plane near the point A can be found from the 
following consideration. As it was established above in the vicinity of the transition line the gap is proportional to the 
deviation from the line. This is valid for any direction of deviation except the direction at a tangent to the transition 
line. Thus, in the vicinity of the point A on the line h x — h xA {A), the gap is 

m~h xA (A)-h xc (h z ,A) (48) 

On the other hand, the gap is given by Ea. l|47|l . Equalizing these two expressions for the gap we obtain the equation 
for the transition line in the vicinity of the point A as 

h xc (h z ,A) = h xA (A)-f(A)h 2 z (49) 

where the function /(A) is generally unknown and can be found numerically only. 

Summarizing all above, we conclude that the MFA correctly describes the critical properties of the transition line 
and determines the transition line with high accuracy. This is because the MFA gives the exact ground state on the 
classical line, which is close to the transition line. Besides, the MFA is asymptotically exact in the vicinity of the 
point F. Therefore, for any value of A, the accuracy in determining of the transition line drops as one moves from the 
point F to the point A. The MFA quality for the case h z = was investigated ini&ii, where it was shown that the 
accuracy of the MFA is high for A > —0.5, and the MFA fails in the limit A — ► — 1 (where the accuracy decreases to 
20%). Besides, the high accuracy of the MFA in the vicinity of the transition line is confirmed by DMRG calculations 
(see Fig|3]and FigQJ. The MFA qualitative correctly describes the line h x = (no gap and sound- like spectrum for 
h z < 1 + A), but one can not expect that the MFA gives correct critical exponents in \ow-h x region. 



IV. THE LOW-li, REGION 

On the line h x = the model Q reduces to the well-known exactly solvable XXZ model in the longitudinal 
magnetic field. In this model three phases exist in different ranges of the magnetic field h z : the ferromagnetic (F) 
phase at h z > 1 + A; the antiferromagnetic (AF) phase at < h z < h z o(A) (h z o(A) is a lower critical field^) and 
the critical phase at < h z < 1 + A (|A| < 1) and h z0 (A) < h z < 1 + A (A > 1) (see FigEJ). 

In the F phase the ground state is saturated ferromagnet M = (S*) = 1/2 with a gap in the spectrum. In this 
region the appearance of the transverse magnetic field does not cause noticeable change in the system properties. It 
results in appearance of a uniform magnetization in the X direction and small decreasing of the magnetization in the 
Z direction. 

In the AF region the system is in a gapful phase with the long-range Neel order M* t and zero uniform magnetization 
M = 0. Due to the gap in the spectrum the effect of the X component of magnetic field in the AF region can be 
obtained in the frame of a regular perturbation theory in h x . The estimate of the first and the second orders in 
h x indicates the appearance of the uniform magnetizations in both the X and the Z directions and the staggered 
magnetization in the X direction as 



(Si) ~ h x + (-i) n h z h x 

(S*) ~ (-l) n M* st + h z hl 



(50) 
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FIG. 5: The ground state phase diagram of the XXZ model in longitudinal magnetic field. 



As follows from the last equations, in the case h z = the applied transverse magnetic field does not cause the 
uniform magnetizations in the Z direction and the staggered magnetization in the X direction^. 

The critical phase is characterized by non-zero magnetization < M < 1/2 in the ground state and by the massless 
spectrum. The low-energy properties in this phase are described by a free massless boson field theory with the 
Hamiltonian 



Ho = \ J dx [(d x Q) 2 + (d x $) 2 ] 



(51) 



where <&(x) and 0(x) are boson and dual field, respectively and v(A,h z ) is the renormalized spin-wave velocity. 
The spin-density operators are represented as2^ 



S z „ ~ M+— !— &.$ + ai(-l) n cos( — +2ttMx 
2itR y ' \R 

SI ~ b Q (-1)" cos (2ttRQ) + &x cos (27i\R6) cos (- + 2ttMx 



(52) 



where ai, b Q and b\ are some constants^ and we identify the site index n with the continuous space variable x. The 
magnetization M(A, h z ) and the compactification radius i?(A, M) are functions of A and h z and can be determined 
by solving Bethe-ansatz integral equations2§i2£. 

Both terms of operator S x in Eq. (|52|l are oscillating when M ^ and are not relevant to the uniform JT-componcnt 
of the magnetic field. But as was shown inSS the second term in Eq. Q52J1 corresponding to perturbation 



Vq = Kbx cos (27ri?6) cos ^ + 2ttMxJ 
has conformal spin S = 1 and generates two other perturbations with zero conformal spin: 



(53) 



Vi = gi cos (47ri?6) 
'2$ 
~R 



V 2 = 02 COS 



AttMx 



(54) 



with gi ~ h 2 , and g 2 
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The scaling dimensions of perturbations Vi, V% are 2r\ and 2/rj (rj — 2ttR 2 ), respectively. The perturbation V2 
describes Umklapp processes and it is responsible for the gap generation in the AF region, where M = 0. But in 
the critical region (h z > 0) the magnetization M ^ and the perturbation V2 as well as the operator Vq does not 
conserve the total momentum and will be frozen out. 

Therefore, in the uniform longitudinal magnetic field h z the critical exponent for the mass gap is determined by 
the only non-oscillating perturbation V\: 

m ~ h]h (55) 

We note that in the special case h z = (|A| < 1) the magnetization M = and all perturbations Vq, Vi and V2 
are non-oscillated 2 «. In this case in the region A > cos[7r-\/2] ss —0.266 the perturbation Vq becomes most relevant 
and determines the mass gap as^ 



m ~ ht"~ 1/r! (56) 

The perturbation V\ corresponds to the spin-nonconserving operator YK^n^n+i ~ Sn&n+i) - 24 This means that the 
behavior of the system at small h x is similar to that of well studied XYZ chain in magnetic field h z with small 
anisotropy in XY planed. Using the results of RefiSi we conclude that in the region A > 1 the presence of small 
X-component of magnetic field leads to a sequence of three transitions with increasing longitudinal magnetic field 
h z . The first one occurs at h z o(A) (77 = 2), where the AF phase transforms to the incommensurate critical (IC) or 
'floating' phase (see FigEJ). In the IC phase the perturbation V\ is irrelevant and the spectrum remains gapless. Here 
the correlation functions display a power-law decay with a magnetization dependent wave vector. There is no LRO in 
the IC phase, except uniform magnetizations M and (S%) = Xxh x (the susceptibility \x is finite in the critical phase). 
Unfortunately, neither the field theory approach nor the MFA allow us to determine the boundary of the IC phase. 
Therefore, this phase boundary is shown on Fig [21 schematically. 

Further increasing of h z leads to the transition of the Kosterlitz-Thouless type taking place at the point h z \{A) 
where r\ = 1. At h z > h z \{A) the perturbation V\ becomes relevant (77 < 1) and the system crosses from the IC phase 
to a strong-coupling regime with the staggered magnetizations in both X and Z directions (AF phase). 

The last transition with further increasing of h z occurs near the point F at h z = h zc (see FigEJ . This transition to 
the PM phase was studied in Sec. III. 

The field theory approach allows also to determine the exponent for the appeared LRO. The transverse magnetic 
field generates the staggered magnetization along the Y axis at |A| < 1 as 

( S n) - ^ ~ (-1)"^ /2 - (-l) n K /(2 - 2V) (57) 

and along X and Z axes at A > 1 and h z i(A) < h z < 1 + A 

(S x n ) - X*h x + (-1)"^ A2 - 2 ") 

(S*) - M + (-i) n fcP-^)/( 2 -»») (58) 

We note, that in the limit h z — > 1 + A (77 — > 1/2), the exponents for the mass gap (|55|l and staggered magnetizations 
in Eqs.Jinj!,© agree with the MFA results in Eas. l{i2l . (|t5j) . lp[j) . 

The staggered magnetization along the Z axis in Ea. H58l) can be derived in the same manner as was derived the 
generated perturbations V\ and V2 in Eas. (|54|) . According to Eas. (|52ll the non-zero contribution to the first order 
correction in h x to the staggered magnetization ((— l) n S*) is given by the following terms in spin-density operators 

(-1)™S* - cos ( ^ + 2itMx ) 



R 



SI - cos {2ttRQ) cos + 2ttMxJ (59) 

Then the leading contribution comes from small distances of order of ultraviolet cut-off (the lattice constant) 

(-l) ai S*(*0 ~ h x [ d 2 z 2e i2 - M ^-^) e xp[-i^il]exp[i^l]cos[2^e(z 2 )] 
J R R 

- h x cQB[2irRQ(z 1 )]~h x (-l) xl S*(z 1 ) (60) 
Thus, the relation between the staggered magnetization along X and Z axes is established 

((-1)"S«> ~ M(-1)»S£> (61) 
which results in the critical exponent for ((— l) n S%) written in Eq.JSSJl. 
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A. Perturbation theory. Mapping to the XY Z model in the magnetic field 



In this subsection we study the behavior of the model in the critical phase from the perturbation theory point of 
view. In some sense this subsection is complementary to the field theory approach. 

For a small X component of a magnetic field the Hamiltonian Q can be written in the form 

H = H + V 

Hq = y^ffn 5*^+1 + S%Sf l+1 + AS Z S Z +1 ) — h z S z 
V = -h x ys x n (62) 



As was said above in the critical region JSJ the model Hq is in the critical regime with non-zero magnetization M. 

The transition operator S x = in V conserves the momentum and changes the total S z by one. Therefore, 

since the unperturbed Hamiltonian Hq conserves the value of S z , the perturbation theory in V contains only even 
orders. 

The operator V acting on the ground state or the lowdying states produces states with 'high' energies E s — Eq = 
e s ^ 1. So, the perturbation theory in V has such a 'twodevel' structure with alternate lowdying and highdying 
intermediate states. As a result the denominator of the 2n-th perturbation order contains n factors with 'high' 
excitation energies e s > 1 and (n — 1) with lowdying excitation e s ~ 1/N. 

Omitting numerical factors, in all orders of perturbation theory one can take out all 'high' energies from the 
denominator and sum up over all highdying intermediate states. As a result, we arrive at the perturbation theory 
with the effective perturbation 

V' = -h 2 x J2KSl (63) 

We note, that the perturbation theory with the perturbation V' coincides with the original perturbation theory 
(|62|l in a sense, that both perturbation series have the same order of divergencies (or power of TV) at each order in 
h x . But numerical factors at each order in h x can be different. 

The behavior of the ground state correlation function { s S x S x +J .) of the model Hq on large distances is known2£ 

<^ +r ) = Mv)^ - Mv) c °l^ M v r) (64) 

and, therefore, due to oscillation of the correlator (S^S^) (for M ^ 0) the sum over n and m can be approximately 
estimated as a one half of the first term: 

E > = J + 2 E (Sf^J ^ ^ + E ( S nSn +1 ) (65) 

The last equation suggests that the perturbation V can be reduced to the operator V" 

V" = -ahlN-bhlY.S x n Sl +l (66) 

with some constants a, b. In order to verify this assumption we compared the matrix elements of the operator V" 
with those of V . We have found that the dependence of the matrix elements of the operators V and V" on N is 
identical. 

As a result of the analysis we arrive at the effective Hamiltonian 

HeS = -aNh 2 x +E((1- bh 2 x )S x n S* n+1 + SlSl +1 + AS*S* +1 ) -^E 5 - ( 67 ) 

Again, the original model l|62|l is equivalent to the effective model ()67|l in a sense, that the perturbation series for 
both models have the same order of divergencies at each order of h z . This fact means that the effective model has 
the same critical exponents as the original model (|62|l . 

One can see, that the effective perturbation V" is nothing but the generated perturbation V\ found in bosonization 
technique Ea. (|54|l with g\ = bh x . So, the above analysis represents another way to obtain an effective perturbation 
responsible for the most divergent part in the perturbation theory. 

Unfortunately, the effective model l|67|l is not exactly solvable one in general case. But it can be used for checking 
the bosonization expression for the mass gap (|55|l in two special cases h z = and A = 0, where the effective model 
reduces to integrable XYZ model and XY model in magnetic field, respectively. In these special cases the exact 
expressions for the mass gap are knowniS^i. They confirm the validity of the critical exponent in Eq. l|55ll . 
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V. SPECIAL CASE A = 1 



When A — > — 1 the point F tends to zero, and exactly on the line A = — 1 both classical and transition lines 
disappear. In the vicinity of the line A = —1 it is convenient to rotate the coordinate system in such a way that the 
Hamiltonian takes the form 

H = H + Vh + VA 

Ho = — s„ ■ s„ + i — h x i)"^ 

V h = -h^S* 

Va = (1 + A)]TS^ +1 (68) 

The unperturbed Hamiltonian Ho describes the critical behavior of the system at < h x < ho, 32 where the estimate 
of ho is ho w 0.53(3)^ The spin-density operators of Ho are related to the bosonic fields by Eas. (|52|) with M = 0, 
where the compactification radius R(h x ) is a function of h x . This dependence is generally unknown, but the limiting 
values of rj(h x ) — 2-KR 2 (h x ) are knowni^ 

V (h x ) = — , h x — > 

T) (h x ) = -, h x ^ h (69) 

The perturbations 14 and Va have scaling dimensions rj/2 and 2r/, respectively. Therefore, according to Eas. l|69|) 
both perturbations Vh and Va are relevant in the critical region < h x < ho and produce a mass gaps 

m h ~ h\l^l^ 

m A ~ ll + AI 17 ^ 2 ") (70) 

In particular, m/, ~ \[h~ z , toa ~ \Zfl + A] at h x — > and to/,. ~ to-a ~ |1 + A| 2/3 at h x —> ho- 

The gap at A = — 1 and for h x ,h z << 1 can be asymptotically exactly described by the spin- wave theory^. The 

validity of the spin-wave approximation is quite natural because the number of magnons forming the ground state is 

small for h x ,h z << 1. The spin- wave result for the gap is 



m= ] Jh z (h z + } ^) (71) 

For h z << h x << 1 Ea. H71[) agrees with the conformal theory result Ea. l|7U|l at r/ — > and gives the preexponential 
factor for the gap. 



VI. DISCUSSION AND CONCLUSION 



As it was mentioned the order-disorder transition induced by the magnetic field has been observed in the quasi- 
one-dimensional antiferromagnet CS2C0CI4 described, as supposed, by the model JD with A = 0.25. It is interesting 
to compare the magnetization curves obtained in neutron-scattering experiment^ with those in the MFA. In this 
experiment the magnetic field H has been applied at an angle j3 ~ 40° to the XY plane. This means that H x ~ H z ~ 
H/y/2. According to£ ^-factors in CS2C0CI4 are g X}V = 2g z = 4.80. Therefore, the ratio of the effective fields in the 
model (QJ is jr — ~ 2. The total magnetization M to t is 

M tot = J^(g x {SZ)+g z (S z n )) (72) 

where (S%) and (S*) are the magnetizations calculated in the MFA at the effective fields h x — 2h z . The AF ordered 
moment is given by M st — g y (j,t,{\(—l) n S%\). On FigEJwe plot M to t and M st as the functions of the magnetic field H 
(H = g h ^ b J )- These magnetization curves are qualitative similar to the experimental ones (Figs. 12 and 14 in Refi£). 
The maximal value of staggered magnetization 1.7/^f, agrees with the experimental magnitude of the AF ordered 
moment 1.6/x;,. The total magnetization at H — > 00 in the MFA M to t = 1-9/ib is consistent with the saturation 
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moment 1.7 Hb estimated in£. At the same time, there is essential difference in the low- field behavior of M st . The 
experimental AF ordered moment is finite at H — 0, while M st — at H = on FigEJ This difference is due to weak 
interchain couplings, which form the magnetically ordered ground state in the real compound at H — and this effect 
is absent in the one-dimensional model £Q) . The behavior of the magnetizations in the vicinity of the critical field on 
Figgis similar to the experimental curve. But the value of the critical field on Figgis 1.6T, while the experimental 
value is 2. IT. We do not believe that the MFA is the reason for this discrepancy. The possible reason may lie in 
the fact that CS2C0CI4 is described by the s = | antiferromagnetic model with strong single- ion anisotropy and its 
reducing to the s — ^ model is approximate. 

Magnetic measurements in£ were performed in the field applied at a fixed angle (3 to the anisotropy axis. It is 
interesting to consider how the properties of the model is changed when the magnetic field is turned from purely 
longitudinal (f3 = 0) to transverse direction ((3 = ^). If the effective field h is in the range (1 + A) < h < h x A(A) for 
|A| < 1 or H x a(A) < h < (1 + A) for A > 1 then there is the critical angle fio, which is defined by the intersection 
of the circle h x + h 2 z — h 2 with the transition line. At this angle the phase transition from the AF phase at (3 > (3$ 
(f3 < flo) for A < 1 (A > 1) to PM phase at (3 < (3q {(3 > (3o) takes place. The staggered magnetization and the gap 
vanish at (3 — (3 as M st ~ |/3 — Po] 1 ^ 8 an d m <~ \(3 — (3 \. 

In conclusion, we have studied the spin-^ XX Z Heisenberg chain in the mixed longitudinal and transverse magnetic 
field. It was shown that the ground state phase diagram on the (h x , h z ) plane contains the AF and the PM phases 
separated by the transition line. The transition line was determined using the proposed special version of the MFA, 
which reduces the XXZ model in the mixed fields to the XY model in the uniform longitudinal field. The MFA 
gives the transition line with high accuracy at least for A > —0.5. This fact is confirmed by comparison of the MFA 
results with DMRG calculations. The MFA gives satisfactory description of the whole phase diagram, though the 
critical exponents of low-h x dependence of the gap and the magnetization can not be found correctly in the MFA. 
These exponents have been found with use of the conformal field method. 

The field theory approach also shows that in the region A > 1 the phase diagram contains IC or floating phase 
characterized by the gapless spectrum and power-low decay of correlation functions. The form of the boundary of 
IC phase can be determined numerically only. But this boundary is located certainly on the left of the classical line, 
where the spectrum is gapped. 

We believe that the modified MFA is suitable for the studying of the magnetic phase transitions of the quasi-one 
dimensional anisotropic magnets induced by the applied magnetic field. The important problem is to take into account 
effects of inter-chain interactions^. 
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